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We discuss an action in which the fermion matrix has single level stout smearing for the hopping 
terms together with unsmeared links for the clover term. With the (tree level) Symanzik improved 
gluon action this constitutes the Stout Link Non-perturbative Clover or SLiNC action. To cancel 
0{a) terms the clover coefficient, Csw has to be tuned. We present here preUminary results of a 
non-perturbative determination of Csw using the Schrodinger functional and as a by-product also 
a determination of the critical hopping parameter, fQ. A determination of the renormalisation 
constant for the local vector current is also given. Comparisons of the results are made with 
lowest order perturbation theory results. 



The XXVI International Symposium on Lattice Field Theory 
July 14-19 July 2008 
Williamsburg, USA 

* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Non-perturbative determination ofcswfor SLiNC fermions 



R. Horsley 



1. 0{a) Improvement 

When constructing a lattice QCD action, even the simplest gluon Lagrangian action has only 
0{a^) corrections. The naive fermion action also has 0{a^) corrections, but suffers from the 'dou- 
bling problem' describing 16 flavours in the continuum limit. A 'cure' is to add the Wilson mass 
term, so 15 flavours decouple in the continuum limit, but the price is that there are now 0{a) 
corrections, so that for example for a ratio of hadron masses 

= ro + ari + 0[a ) . 

THh' 

The Symanzik approach is a systematic improvement to 0{a") (where in practice n = 2) by adding 
a basis (an asymptotic series) of irrelevant operators and tuning their coefficients to remove com- 
pletely 0(a"^^) effects. Restricting improvement to on-shell quantities the equations of motion 
reduce the set of operators in both the action and in matrix elements. Indeed, for 0{a) improve- 
ment only one additional operator in the action is required 

■^clover 

the so-called 'clover term'. So if we can improve one on-shell quantity this then fixes Cgw as a 
function of the lattice spacing a or equivalently of gQ, so that all other physical on-shell quantities 
are automatically improved to 0(a), i.e., we now have 

^ = ro + 0(a^). 
niH' 

Matrix elements still require additional 0{a) operators, for example 

<<4 = {\ + bAamq){A^+CAad'jfP) 
gfi = {l+bpamq)P, 

with 

An easily determined quantity is the quark mass, determined from the PCAC relation^ 

wj {d{f{Ao{xo) + CAad,fP{xo))0) 
" 2{P{xo)0) 

Choosing different boundary conditions or operators, O, gives different determinations of the quark 
mass mq'^'\ i = 1, 2. If the quark mass is improved then its errors are O(a^). So we can determine 
improvement coefficients, Csw, • • by finding the point where 

w;(l) wi(2) 
ruq^ ' = mq"^ ' . 

^This is equivalent to considering the renormalised quark mass = rrH^' as the difference is just a 

numerical factor, which in the chiral limit does not effect considerations of O (a) -improvement. 
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The ALPHA Collaboration achieved this by means of the 'Schrodinger functional', [|l|]. Dirichlet 
boundary conditions are applied on the time boundaries to the fields. For the gluon fields fixing 
them on ;co = and T is then equivalent to a constant chromo-electric background field (which 
means that simulations with niq ^ with no zero mode problems are possible), while the fixed 
quark fields (p,p) can be taken as sinks/sources to build operators for correlation functions. For 
example here we can take at the lower boundary xo = 0{i = \) and upper boundary xo = T {i = 2) 




So we can investigate PCAC behaviour at different distances from the boundaries. Redefining the 
quark mass slightly (but in a way which coincides to 0{a^) in the improved theory, |Q]) to eliminate 
the unknown ca (m^' —^ M) we can define improvement when 

(M,AM) = (0,0), 

where 

are chosen at some suitable ;to, [@]. This gives the required critical c*^ and K*. 

There are (small) ambiguities due to the finite volume used. In an infinite volume we expect 
O(aAgco) contributions (in the chiral limit, otherwise there are also extra 0{amq) terms) due to the 
different boundary conditions or operators chosen. In a finite volume there are additional 0{a/Ls) 
terms. 0{aAfx:D) ^ as a (or g^) 0, but 0{a/Ls) ~ 0{\/Ns) (where = aNs). We can either 
keep Ls fixed in physical units as a — > (the 'constant physics condition') so 0{a/Ls) — > 0, or 
alternatively simulate for several values of A^^ and extrapolate to A^^ oo. The 'Poor man's solution' 
is to evaluate at large p —> oo (small a) and subtract this result. Practically we have found that for 
Csw this 0{\/Ns) term is negligible, while for Zy, this subtraction is about a 1% affect. 

2. The SLiNC action 

We shall apply the Schrodinger functional formalism to 2 + 1 flavour stout link clover fermions - 
SLiNC fermions (Stout Link Non-perturbative Clover). In a little more detail 

= L { ^'W{^)U^ (x + /i) [7^ - 1] i/^(x - A) - K-v7(x)f7^ (x - A) [7, + 1] i/a(x + A) 

X 

+"V7(x)v/-(x) + \cs„{gl)'W{x)o^yF^y{x)\ir{x)] . 

The hopping terms (Dirac kinetic term and Wilson mass term) use a once iterated stout smeared 
link or 'fat link', 

tj^ = exp{/e^(x)}[/^(x) 
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(V^ is the sum of all staples around U^) while the clover term remains built from 'thin' links - they 
are already of length 4a and we want to avoid the fermion matrix becoming too extended. Smearing 
is thought to help at present lattice spacings and the stout variation is analytic which means that 
the derivative can be taken (so the HMC force is well defined) and perturbative expansions are also 
possible, 1^. 

To complete the action we also use the Symanzik tree-level gluon action 



Sg 

together with 
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3. The lattice simulation 



The lattice simulation used the Chroma software library, [Q]. The Schrodinger Functional 
details follow |j5|]. All results were generated on 8^ x 16 lattices using the HMC algorithm. A 
mild smearing of a = 0.1 was used. A series of simulations were performed (typically generating 
0(3000) trajectories), quadratic and then linear interpolations of the M, AM results being used to 
locate the critical point. 

We thus have a two-parameter interpolation in c^w and k: which is split here into two separate 
interpolations. First plotting AM against M and then interpolating to M = gives AM {csw, 'Q(c,vw))- 
A typical result is shown in Fig. |l|. These values of AM for M = for various j8 values are then 
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Figure 1: AM against M for j3 = 6.00 (filled symbols) together with quadratic interpolations to M = (the 
open symbols). 

plotted against c^w as shown in Fig. ^ AM = then gives c*„. 

A similar procedure yields K* : plotting M against 1/k and interpolating to M = gives 
Kc{csw)- Then subsequently plotting AM against I/jq. and interpolating to AM = gives k:*. 
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Figure 2: IsM at M = against c^,,, for various values of /3 (filled circles) together with linear interpolations 
to AM — (open circles). 



4. Results 

The results for c*^, and K* against are plotted in Figs. ^, ^respectively in the range j3 < 5. 10. 
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Figure 3: c*j,, against for various values of /3 (circles), together with a polynomial interpolation (line). 
Also shown is the perturbative result. 

The lowest order perturbative limit has been computed for both c*„. and K*, and is also shown 
in the figures. An interpolation between the numerically determined points is also shown. For 
both c*„, and K* a 6th order polynomial in proved sufficient. (These interpolation functions aie 
constrained to reproduce the perturbative results, [Q], in the j8 ^ oo limit. Therefore, they have four 
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Figure 4: k* against for various values of j3 (circles), together with a polynomial interpolation (line). 
Also shown is the perturbative result. 

free parameters.) This smooth fit between the points gives an estimate for cl^^, which will be used 
in the action for future generation of configurations. 

For c*n, the polynomial only tracks the perturbative solution for small values of gQ. This is 
perhaps not surprising as the tadpole improved, TI, estimate is c^J^, = u[^^ /uq, [^], which is to be 
compared with the unsmeared case of cj^^, = 1 /mq where uq is the average plaquette value and Uq'^ 
is the smeared value. As smearing increases the plaquette value this indicates that c*^ can be large. 
For K* on the other hand as jc" = 1/ (8mq'^^) we expect that it is 1 /8. This is true for reasonably 
fine lattices, however k* does begin to decrease for larger values of ^q. For nj = 2 the same 
phenomenon occurs: for larger g^, K* begins to decrease (after initially increasing). 

Finally in Fig. § we show the vector renormalisation constant. This is computed using the 
vector current in the ratio of a three-point to two-point function where the sinks/sources are built 
using O(') as described earlier in section |I|. 

5. Conclusions 

Non-perturbative 0{a) improvement is a viable procedure for (stout) smeared actions with 
typical clover results being obtained. (Other recent results for 2 + 1 flavours are given in [§, 0].) 
As a decreases we need a significant Csw ^ c'sw = 1 for 0{a) improvement. We are now seeking a 
region where a ~ 0.05 — 0.1 fm. Improvement, which is presumably an asymptotic series, brings 
an advantage for smaller a say a < 0.1 fm. The two extremes for a are simulations at small a with 
'large' nips when there is no continuum extrapolation but a chiral extrapolation, or alternatively 
simulations at 'coarse' a with mp^ ~ rrijc when there is no chiral extrapolation but a continuum 
extrapolation. Of course the Schrodinger functional does not tell us a; for this conventional HMC 
simulations are required. Some preliminary results indicate that around j8 < 6.0 we have a < 
0.07 fm. Final results, including larger lattice size comparisons will be published elsewhere, 
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Figure 5: Zy against for various values of j3 (circles), together with a polynomial interpolation (line). 
Also shown is the perturbative result. 
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